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f~^ I ing the relative acceleration between particles on different space-time trajectories. 
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1 Introduction 

According to the equivalence principle, the motion of structureless test parti- 
cles in a gravitational background field is determined only by the space-time 
geometry: the worldlines are geodesies |jl|. However, as only relative accel- 
erations (tidal forces) have an invariant meaning in general relativity, the 
relevant information about the space-time geometry as a gravitational infiu- 
ence manifests itself in the rate at which worldlines diverge (or converge). 
The equation describing the rate of geodesic deviation in a quantitative way 
is well-known since at least half a century [0], and numerous applications 
have been developed up to the present time !§]-[§]• 

Things become more complicated for test particles which are not struc- 
tureless, but carry e.g. non-vanishing charge or spin. In such cases the world- 
line of the test particle in general is no longer a geodesic, but is modified by 
electro-magnetic and/or spin-orbit forces [0, |^, §]. In these cases one also can 
write an equation for world-line deviations, describing relative accelerations, 
which are now due to the combined effect of gravitational and/or electro- 
magnetic forces. Even a generalized string-theoretical version has recently 
been proposed |]10 . 



Nevertheless, in the presence of such additional forces the geometrical 
interpretation of the world-line deviation is not completely lost. For exam- 



ple, it is well known [1^, |1J], that the geodesic lines in the five-dimensional 
Kaluza-Klein theory represent the worldlines of charged massive particles, 
moving under the infiuence of the Lorentz force in a curved four-dimensional 
space-time. In fact, the equivalence principle holds in five dimensions as 
well, although in four dimensions it applies separately to classes of particles 
characterized by the same value of the charge-to-mass ratio q/m. Recently 



it has been checked [jT5| that also the geodesic deviation equation in the five- 
dimensional Kaluza-Klein manifold leads to a generalized deviation equation 
for the world lines of charged particles, which can be obtained by a direct 
variation of the world line's equation in the four-dimensional space-time. 
Similarly, the world lines of particles with spin can be obtained from the 
supersymmetric extension of the geodesic equation for simple point particles 
[0, 0. In spite of these developments, examples of a thorough analysis of 
solutions to the world-line deviation equations, as well as discussions of their 
physical implications, are scarce in the literature (||16||). 



In this article we initiate a self-contained discussion of the world-line de- 
viation equations in Einstein-Maxwell theory. We first derive the basic form 
of the equation, taking the equation of motion itself as a starting point. 
We also show that these equations can be derived from an action principle. 
We then turn to discuss two examples for which the nature of the solutions 
to the world-line deviation equation can be explicitly exhibited. The first 
case describes the effects of combined non- stationary gravitational and elec- 
tromagnetic fields, in the form of a gravitational pp-wave. The second case 
concerns the analysis of deviation of world lines in stationary spherically sym- 
metric field generated by an electrically charged mass (Reissner-Nordstr0m 
solution). 

In part, the choice of these particular solutions is determined by their 
high symmetry properties, which make it possible to arrive at exact and ex- 
plicit solutions for the background field, for the world lines, as well as for the 
deviation itself. But another reason for investigating this type of background 
solutions derives from the physical aspects: according to General Relativity 
the electromagnetic field acts on massive charged particles in two ways, di- 
rectly via the Lorentz force it exerts, and indirectly, via its contribution to 
the ambient space-time curvature. The combination of these effects may lead 
to interesting new states of motion. Indeed, the Reissner-Nordstr0m solution 
provides a particularly illustrative example of this situation, as the repulsive 
electric force on a test particle can be compensated by the attractive gravi- 
tational force. 

The supersymmetric extension of the models discussed here describes the 
modifications in the physics due to the presence of spin. These aspects will 



be discussed in a separate paper [17 



2 World-line deviation equations 

Before considering the more general case, let us recall the derivation of the 
equation for geodesic deviation relevant to simple chargeless point particles. 
Consider a continuous smooth vector field u'^lx) with the property 

u^'V^u^ = 0. (1) 



A vector field of this kind, being transported parallel to itself, represents the 
field of vectors tangent to a smooth set of geodesic curves x'^{t): 

■■ ^"' (2) 



dr 



) 



where as usual for time-like geodesies we take the parameter r to represent 
proper time. In terms of proper-time derivatives eq.(|I]) then becomes 

D^x^' _ cPx'' ^ dx^ dx" _ 

Now consider a slicing of space-time into smooth (l+l)-dimensional space- 
time planes formed by time-like geodesies, and in each of these planes a 
second smooth set of geodesic curves x^[\] r], connecting the points on the 
geodesies corresponding to a fixed value of the proper-time parameteiQ r. 
Thus, for each value of r there is a smooth curve in the plane parametrized 
by the real affine parameter A. Let n^{x) = dx^/dX be the field of tangent 
vectors to these curves. We observe that 

Combining eqs.(0) and (^) we then obtain 

= n^Vx (u'V^u'') + u^-" [Va, V,] m'^ = R^^fu\^n\ (5) 

In terms of derivations in parameter space (proper time derivatives) this 
equations reads: 

-^^ = R>.Ju\^n^, (6) 

which is the standard text-book result, as can be found e.g. in |^ PT|] , 
stating that the rate of change of the velocity at which geodesies diverge, as 
measured along the curves x^{\) at fixed r, is proportional to the geodesic 
separation, with the proportionality factor given by the components of the 



^In order for these curves to be well-defined, one fixes the origin of proper time on each 
geodesic to be located on a smooth, differentiable curve cutting all the geodesies in the 
slice. 



Rieinann curvature tensor in the direction of the geodesies and the curves of 
equal proper time. 

It is now straightforward to extend this derivation to curves representing 
the world line of charged particles. We start from the equation of motion for 
a particle of mass m and charge q: 

^ = u'^V.u^ = ^F^y. (7) 

The right-hand side of this equation represents the electro-magnetic and 
Lorentz force acting on the particle. Next we obtain an equation for the 
relative acceleration of two world lines, as measured by the change in tangent 
vectors n^ = dx^/dX to curves connecting points of equal proper time on 
smooth planes of world lines. The equation, derived by a series of steps 
similar to those of (^, reads: 

^ = Rx.fu\^n^ + - n^ + - V^n n^n\ (8) 

Ut^ m Ut m 

In practice it is often simpler to work with the non-explicitly covariant equa- 
tion 

n'^ + (2TJu^ - ^F^)i if + (u''u''d,T,^ - ^u^O.fA n" = 0, (9) 
V m / V 771 J 

which is obtained from (H) by substitution of the explicit expressions for the 
covariant derivatives and the Riemann curvature tensor. 

We observe, that this equation only determines the rate of divergence of 
world lines of particles with the same charge-to-mass ratio q/m. For neutral 
particles equation (|^) remains valid as the special case for g = 0. Eq.(|) does 
not apply to the relative acceleration between particles of different charge- 
to-mass ratio, such as a neutral and a charged particle. On the other hand, 
it is clear that a neutral particle can feel the presence of an electro-magnetic 



field through the influence of its accompanying gravitational field |jTT . 

We also observe, that the direct influence of the electro-magnetic field 
on charged particles is strictly linear in the field strength F^,^, with coupling 
constant q/m. No terms quadratic in F^j, are present in this equation in the 
context of classical general relativity. 



Observe, that eqs.(|^ and (H) can be combined to show that u ■ Dn/Dr 
is a constant on geodesies: 

— [u =0. 10) 

dT \ DtJ ^ ' 

Actually, this equation can be integrated to give the stronger condition 

Dn , . 

«.- = 0. (11) 

This follows, if we consider two worldlines Xi{t) = x^{t;X) and X2{r) = 
x^(r; A + AA) ^ x1{r) + As^(r), with 

Ax'^(r) = n^(r)AA. (12) 

As both X(^i^2){t) are solutions of the equations of motion, it follows that 

gfiuixi)u1u'( = -1 = g^u{x2)ui^U2 

~ g^,u{xi)u1u'( + 2AXg^^{xi)u'^--—. 

JJt 

As a result only three components of the rate of change of the deviation 
vector n^{T) are independent. 

The equation (H) for the world-line deviation of charged particles in gen- 
eral relativity, and per force the geodesic deviation equation (|^), can be 
derived from a covariant action principle. The relevant lagrangean reads 

2 ^^ Dt Dt 2 ^ 2m ^ Dr 2m ^ 

(14) 

After some manipulations involving the Bianchi identity for the electro- 
magnetic field, the world-line deviation equation (^ is now obtained by 
requiring the action S = J Ldr to be stationary under variations of the 
deviation vector n^: 

^ = 0. (15) 

This action can serve as a starting point for a Hamiltonian formulation and 
for quantization. 



3 Motion in PP-waves 

3.1 Particle worldlines 

The combined Einstein-Maxwell equations admit parallel plane-wave solu- 
tions, known as pp-wave^. The non-singular metrics of these solutions are 
represented by the line-element 

c^dr'^ = —dudv — K{u, x, y) du^ + dx^ + dy^, (16) 

which reduces to flat Minkowski space-time when the metric Mw-component 
vanishes: K{u,x,y) = 0. As is implicit in (|16]), we choose to work with the 
light-cone co-ordinates 

u = ct — z, V = ct + z. (17) 

As the determinant of this metric is negative definite, it is manifestly invert- 
ible everywhere. 

The solutions for the vector potential of the electro- magnetic field are of 
the form 

A, = 6^^AM + S^A^iu). (18) 

The corresponding electric and magnetic fields are wave-like, propagating in 
the 2;-direction, transverse and orthogonal: 

— = -eijBj = A[{u), i,j = {x,y). (19) 

Here the prime denotes a derivative with respect to u. These fields satisfy 



the free Maxwell equations in a space-time with metric (16). The Einstein 



equations, with the energy-momentum tensor of the electro-magnetic field, 



reduce now to the single p.d.e. [19, O, K 



A^ransK = ^-^ El (20) 

Here Atrans is the ordinary 2-dimensional laplacian in the x-y-plane. The 
general solution for the metric component K is 



AneoG 



K{x, y, u) = -^ Ef {x' + y') + /(«, C) + /(«, C), (21) 



^4 



•^Some reviews can be found in |13, Oq, ^\ [T2|. 



with /, / conjugate holomorphic functions of the complex transverse co- 
ordinates 

( = x + iy, ( = x-iy. (22) 

However, as constant or hnear functions of {(, () do not give rise to a true 
gravitational field (the full Riemann tensor vanishes for them), we take f{u, () 
to be an analytic function of the type 

oo 

/(n,C) = E'^nWC" (23) 

n=2 

Because of its behaviour under rotations in the transverse plane we can iden- 
tify the first term with n = 2 as free quadrupole gravitational waves. In 
contrast, the first term in (PT|), which is the special solution due to the 
presence of the electro-magnetic waves, is invariant under rotations in the 
transverse plane. 

The equations of motion (|^, adapted to this case, can be obtained from 
the world-line lagrangean 

1 n 

L = - (-iiv - Kii^ + x'' + f) - - {A^x + Ayv) . (24) 

Three of the equations of motion can be written as 
ii = constant = 7, 



X = - — K^^ Ex, y = - — K^y Ey. 

2 mc 2 mc 



(25) 



The fourth equation, for the light-cone co-ordinate v, can be replaced by the 
conservation of the world-line hamiltonian 

(26) 







H = - 


2p,pu + 2Kpl + - { 


J^l 


+ na, 


where the 


(covariant) momenta are defined by 








Pv 


1 
= -2" 


7 


1 


- Kii, 




n. 


= Px + 


q 

— A^=X, Uy = Py + 


m 



(27) 



Here {px,Py) are the standard canonical momenta in the transverse plane, as 
used to define the Poisson brackets. It is convenient to use the first equation 
( P^D to change the world-line parameter from proper time r to light-cone 
time u. Furthermore, equation (^) can be solved for the velocity component 
i), or equivalently v' = 7?): 



2H 

-2 



v' = {x[f ~ K - —, (28) 



the last term being a constant of integration. In terms of the constants 
H and p^ the dynamical problem is now reduced to that of solving for the 
motion in the transverse plane. Consider the case of electro-magnetic and 
pure quadrupole free waves; writing /€+ = 4Refi;2, Hx = —4:ImK2 the metric 
component is 



— T- ^i (") (^ +y ) +^~ [X ~y ) + Kx iu) xy = - 



This metric allows additional Killing vectors in the x-y-plane, making the 
dynamical problem fully solvable in terms of first integrals of motion. The 
additional constants of motion are given by 

J = J^Ux + JyUy + 2p^ ( j;.x + J^y) + A 

(30) 

= J^x - .hx + Jyij - jyy + A, 



where the co-efficient functions J(x,y) (u) and A are the solutions of the dif- 

dj 
dr 



ferential equations resulting from the requirement ^^ = 0: 



j; + l{-. + '-^Ef)j, + ^j.. = 0, (31) 

A' = (Ex-Jx + Eyjy) . 

mc 
These equations can be summarized in the compact form 



They admit a two-parameter set of solutions, which imphes the existence of 
two hnearly independent constants of motion fixing the values of the velocity 
components x and y. 

This can be inferred from the trivial vanishing of the Wronskian W{f, g) = 
if 9 ~ Iq') of ciny two solutions, J(i) and J(2), of (|30D , which reduces the orig- 



inal set of four constants of integration of the system (^) to two constants 
only. 

3.2 World-line deviation 

We now discuss the generalized world-line deviation equation (H) in the con- 
text of the electro-magnetic and free quadrupole gravitational waves de- 
scribed by the metric ([T6|) , (^) and the electro-magnetic vector potential 

The only non- vanishing components of the Riemann curvature tensor in 
this case are 



__1 _ _\ ( STreoG ^ 

D ^ Ty- ^ f StIEqG 2\ foo\ 

Ryuyu = -^Kyy = -^ (^~'^+ "^ ^^ ^i ) ' ^'^'^^ 

whilst the non-vanishing electro-magnetic field strength has two non-zero 
component 

^ ux ~ ^X) uy ~ y V '^/ 

It is then straightforward to write down the explicit form of the world-line 
deviation equations. For the u component of the deviation vector it reduces 
to 

^ = ft" = 0- (35) 

Hence the rate of change of n" is constant: n^ = /3 = constant, and 

n'^ir) = n"(0) + (3t. (36) 



The world-line deviation in the transverse x-y-plane is described by the equa- 
tions 



D^n' ..i 1 dKi ^ ^^ .„ 

n + - —r- ^ + 7 ^,j^ 



Dt^ 2 dr 

„2 

A ,. , .X _ 

2 '^ 2 '^^' ■" mc '■" mc 

This can be simplified, using u as independend variable, to read 



(37) 



- — Kim^ + -Ki.x^n'' - — E'n" - —Eih''. 



u\' 



(38) 
This is a system of coupled ordinary linear second order differential equations 
with field dependend co-efficients. In spite of their complicated form, the 
equations for the transverse components n^ and n^ can be solved in terms of 
the constants of motion (|30D, (31). First define the linear combinations 



Ni = ni- x'y + Xi(n")', i = (x, y). (39) 

By substitution of the explicit form of the equations ( pH] ) we obtain 

N: = -Ik,,N,. (40) 

We observe, that these equations are identical with eqs.(|3l|), (|32|) . Therefore 
they admit the solutions 

N^ = aJ^, Ny = aJy, (41) 

where a is an arbitrary constant of proportionality; this constant serves to 
normalize the solutions for n^, and can be absorbed into the parameter A 
measuring the deviation between the world lines. Together with the solution 
(|36D for n" we now have a full solution for three of the four components n^. 
The last remaining component is the second light-cone component n^, 



10 



which is subject to the equation 
1 D^n'' 



72 Dt^ 



[rf)' + K'n'' + Kiv}\ + iT'fn")' + Kin'' + - K^n'' 



Kijx"x''n'' - Kijx^'n^ 



2q_ 

mcy 



Ei ( n" + -KijX^n'' 



2q_ 

mcy 



E'x"n''. 



(42) 

Upon substitution of the definitions (^) once more, using the deviation 
equations (|35D and (^0]), and using the equations of motion (|25|) to ehminate 
the exphcit dependence on the electromagnetic fields, this equation can be 
cast into the form 



{n^y - 2 {N[xii +{K- (xO^)' n- - {k - {:^^ 



'^^^n^' 



0. 



(43) 



Finally equation (|28|) can be substituted, which together with equation (^) 
gives the final result 

[n" + 2t;(n")' - t;V - 2N^x,]" = 0. (44) 

It follows that the second light-cone component of the world-line deviation 
is given by 

n'' + 2t;(n")' - f'n" = 2aJ-Xi + A + Bu, (45) 

with A and B constants of integration to be fixed by the initial conditions. 
Observe, that starting out with n" = (n")' = 0, the right-hand side directly 
gives n" itself. 



It is well known |2y] that the pp-wave solution may be encoded in another 
frequently used co-ordinate system, in which the line element takes on the 
following form: 

ds'^ = dudv + g22{u){dx^Y + g^-i{u){dx^Y , (46) 

where the light-cone variables u and v are defined as before. 

The space-time endowed with metric (11) admits Petrov's G^ group of 



isometrics |T8| , |T9|. The three Killing vectors, generating the Abelian sub- 
group of G5 are defined as: 






(47) 



11 



The first Killing vector is covariantly constant, isotropic and orthogonal to 



second and third vectors in (47). The Lie derivative of the metric along any 
Killing vector is equal to zero: L^{g^^) = 0. 

Using these symmetries, we can explicitly integrate the world-line equa- 
tions, and then also obtain the explicit solutions for the deviation 4- vector in 
these co-ordinates. The details of these calculations can be found in ([0). 

4 Reissner-Nordstr0in fields 

4.1 Worldlines and integrals of motion 

The Reissner-Nordstr0m solution of the Einstein-Maxwell system of equa- 
tions describes the external gravitational and electro-magnetic fields of a 
static and spherically symmetric distribution of mass and charge, or — in the 
case one takes it as a complete solution — a static charged black hole of mass 
M and charge Q. The gravitational field is described by the line-element (in 
natural co-ordinates with c = 1) 



g.Jx^'dx" = -dr^ = -B(r)dt^ + —— dr^ + r^ idO^ + sin^ 6 dcp^) , (4J 

B{r) ^ ^ 

with 



2M Q^_f-. M\^ Q^-M'' 
The corresponding solution of Maxwell's equations in this space-time is 

- — dt, F = dA= - — 
Airr Anr^^ 

Details of this solution of the field equations are discussed in standard text 
books, see e.g. [^. In the following we assume M^ > Q^ . 



BW = 1- — + ^= 1-- W^^^^. (49) 



A = Aadx'' = — — dt, F = dA= —^ dr A dt. (50) 



We immediately proceed to the solution of the world-line equations for a 
test particle of mass m and charge q. As the spherical symmetry guarantees 
conservation of angular momentum, the particle orbits are always confined to 
an equatorial plane, which we choose to be the plane 6 = 7i/2. The angular 
momentum J is then directed along the z-axis. Denoting its magnitude per 
unit of mass by £ = J/m, we have 

dip 

dr r 

12 



(51) 



dt 


qQ 


^ A-Kmr 


dr 


~ I 2M : Q2 

r r^ 




In addition, as the metric is static outside the horizon r+ = M+y/M'^ — Q'^, it 
allows a time-like Killing vector which guarantees the existence of a conserved 
world-line energy (per unit of mass m) e, such that 

(52) 

Finally, the equation for the radial co-ordinate r can be integrated owing to 
the conservation of the world-line hamiltonian, i.e. the conservation of the 
absolute four-velocity: 

A-Kmr J \ r ^'^ j \ "^^ j 

lYiovn this we derive a simplified expression for the radial acceleration: 

dr'^ ATimr^ \ A-Kmr J r^ \ ' r" \ r r"^ j 

(54) 
Eq.(]53|) can in principle be integrated directly. However, to get directly an 
approximate solution to the equations of motion one can also study per- 
turbations of special simple orbits. We follow both paths and compare the 
results. First, we find a solution to the equation for bound orbits in terms of 
an integral, from which we can compute e.g. the periastron shift, similar to 
the case of the Schwarzschild solution; secondly, we can solve the simple case 
of circular orbits, and then study the problem of generic bound orbits by 
considering the world-line deviation equations for the special case of world 
lines close to circular orbits. This is the subject of the final section of this 
paper. 

For generic orbits, we first construct the orbital equation from the sepa- 



rate equations of motions for r{T) and (/'(t), eqs. (|5lD and (^ 

{AKmM - eqQ) 




e^ 




2Kmr 



(55) 



\AKmJ I r^ 
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We parametrize the solutions of this equation for £ 7^ as quasi-Kepler orbits 

^^^^ = l + eZsyi^y ^^'^ 

with e < 1 for bound states. The function y{<f), which is hnear for Kepler 
orbits, here remains to be determined. As the extrema of the orbit (the apo- 
and periastron) are reached for y{ip) = (2n + 1)71 (apastron) and y{ip) = 
27in (periastron) eq. (|55|) evaluated at these extrema give two independent 
equations relating e and tq to the other parameters: 



e^-1 -- 


{AtiitlM - eqQ) 
27fmrQ 


k^'^'^'- 


( ^^ y] 


(l + e^ 




2Mf / ^ 2\ 




+ e*), 




and 










(47rmM - 

2711717 


'^o V 








4 ^ ' 







(57) 



(58) 



With the help of these relations we find for the function y[}p) the first-order 
differential equation 



' dyV 
dLp) 



1 + ^2 


[i-(/)1 

V47rm/ 


6M Q2 

ro r^ 


2e / ^ 2gn e2Q2 

M cos y H 5— cos y 

ro V ro ; rg 


A + Sco 


5 y + C cos^ y 





6 + e^ 



(59) 



It follows that the total change in the orbital angle </? between two periastrons 
is given by 

Av^ = /" ^" (60) 

h V^ + Ecosy + Ccos^y ^ ^ 

14 



As a result we obtain an expression for the periastron shift per one revolution: 
., = A,-2..24--^Uof^^^.S;V (61) 



4.2 World-line deviation near circular orbits 

We observe that for circular orbits r = R = constant, the expressions for 
dr/dr, eg. (p^D , and d'^r/dr'^, eg. (|5^) , must both vanish at all times. This pro- 
duces two relations between the three dynamical guantities {R, e, i), showing 
that the circular orbits are characterized completely by specifying either the 
radial co-ordinate, or the energy, or the angular momentum of the test par- 
ticle. In particular, the eguation for constant radial velocity gives: 

qQ \' ( 2M Q'\( f\ 

Substitution of this result into the expression for the radial acceleration, 
multiplying by i?^, and eguating the result to zero, leads to : 

e _ 

R 




(63) 



\4:7rmJ \ R^ \ R R^ 



Note that the right-hand side vanishes for chargeless test particles, i.e. on 
geodesies. In the limiting case Q = we reproduce the well-known result 



MR'-e{R-?,M) = Q ^ i?= _ |l + ^/l_l^), (64) 



leading to the reguirement R > 6M for stable circular orbits to exist pO . 

Because of the spherical symmetry and the conservation of angular mo- 
mentum all test-particle orbits are planar. Therefore the deviation of any 
bound orbit from a circular one can always be computed with reference to a 
circular orbit in the same plane, which may be chosen to be the eguatorial 
plane 6 = 7r/2. For such world-line deviations the component of the deviation 
out of the plane is always zero: n^(r) = 0. 

15 



On the other hand, one might also be interested in the motion of a cloud of 
testparticles, as for example in disks of matter surrounding stars or planets, 
whose orbits are close but not necessarily precisely in the same plane. In 
such a case the deviation between the planes of the orbits, parametrized by 
n^ , oscillates between positive and negative values. 

In order to analyse the deviations quantitatively, we use the connection 
coefficients and curvature components as listed in the appendix, together 
with the expression (|50|) for the electro-magnetic field strength, to study the 
world- line deviation as given by eq.®, with reference to a nearby circular 
orbit of radius r = R. We begin with the equation for the component n^ out 
of the equatorial plane. Using the properties of the circular reference orbit: 

r = R = constant, u'" = 0, 

(65) 

whilst F^j^ = 0, one finds from the deviation equation in the form (P) 



n 



+ {u^fdeT'n' = h' + u^n' = 0, (66) 



with 

'^^ = {u-f = ^. (67) 

We see that indeed either n^ = 0, or else n^ oscillates at frequency / = uj/27r, 
as the orbit is tilted with respect to the equatorial plane. 

Next we solve for the deviations from the circular orbit within the plane. 
For the components (n*, n"^, n"^) the deviation equations take the general form 

n* + Yjh^ + m'^jn^ = 0, i,j = {t,r,ip). (68) 

The coefficients in this equations have the structure 

/O7*, 0\ /000\ 

7 = 7^t 7; , m = m, . (69) 

Vo7^o; vooo/ 

The precise form of the matrix elements follows from eq.(|D or (P); they will 
be given shortly. Now our treatment applies to orbits deviating little from a 

16 



circular orbit, which still decribe bound states. Therefore r = R + Ar, with 
Ar = n'^AX, must remain bounded, e.g. we look for oscillating solutions of 
(pSI). From eqs. (|55D, ( |5U[ ) it then follows that with standard initial conditions 
(clocks and angles synchronized in the initial apastron: n*(0) = n'^(O) = 
0,n^(0) = maximal), the solutions take the form 



n* = nQsincJiT, n^' = UqCOSujiT, n''^ = tIq sin ujiT. 



(70) 



With this Ansatz, eq. (|68D can be written more explicitly as 

T)ni = 0, 
with the matrix T defined by 



(71) 



/ 



-UJ, 



-Ui 



(2r,, 



■m r 










V 










np 



-UJ^ 



J 



(72) 



The solutions for the frequency ul follow by solving the characteristic equa- 
tion det T = 0: 



LoUC + ^a + ul 



with the coefficients defined by 

c = (-9.r,/ + 4r,/r,/) (« 



m 



t\2 



0, 



(73) 



"^rJ- ipip ^ ^^ rip ^ ipip 



p\2 



(U^) 



m m 



\ 1 / ryy.^ 



(74) 



Disregarding the trivial solution cu^ = 0, which occurs with double degener- 
acy, the solution of (|73D of interest to us is 



uo. 



C+—a 
m 



(75) 
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We compute the values of a to zeroth, and C to first order in qjm. The 
connection coefficients F^- ^ have been collected in appendix (0.a). The elec- 
tromagnetic field strength is given in eg. (|5D[) . Finally, eqs. (^) and 
toegether with (|6^), imply for the velocities 



R^^e _, e (^ , qQ 



where we have neglected terms of higher order in the small parameters 
{M/R,Q/M,q/m)] and 

K)' = ^' = ^- (77) 

As a result of the various substitutions one obtains 



C = -u 



■^ _ 6M _Q^ _ 2qQ 
R MR inrriM 



(78) 



q 2qQ ^ 

— a = —- 7T^ + 

m AmnM 



Observe that all terms of order qQ/AiimM cancel. Indeed, deviations in the 
orbital period are not expected from the Coulomb interaction, but only from 
relativistic corrections to the Coulomb interaction, which are of higher order 
in {M/R,Q/M,q/m). 

The solution of the characteristic equation in this approximation thus 
becomes by 

which gives for the frequency 

/ 3M Q^ \ 



To complete the solutions ([70D, eq.(|TI]) establishes a relation between the 
components n^^, which for the solutions (^) becomes 



This guarantees that the perturbed orbit is again a solution of the equations 
of motion. We can interpret this in terms of observable co-ordinate differences 
by writing Ax^ = n^/S.\. 

We finish with some physical observations relating to our solution ([fO|) 
with frequency (|80[) . First, note that there are different solutions uji for 
different values of the test charge q: there is the frequency ojq for neutral test 
particles (g = 0), a larger one u;+ for test and central charge of equal sign, 
and a smaller one uo^ for test and central charge of opposite sign. 

If we compare our approximate solution with the parametrized exact so- 
lution (|6]): 

r(v9) ^ i? + Ar ^ i? - eRcosy{ip), (82) 

we find, using the initial conditions in the apastron {y = —it), that to first 

order 

(Ar)o n^AA 

Thus our solution describes an approximate ellipse (Kepler orbit), with ec- 
centricity e. 

We can also calculate the precession of the periastron. Observe that as 
the orbit reaches its extremal radius, the shifts in orbital angle and time 
vanish: n^^^j, = ±nQ at UiT = kii, with k integer; therefore nl^^^ = n^xtr = 0- 
Now we can calculate the period between two periastra of the perturbed 
orbit: 

T = dr^ = / dr u' + h'AX". 

Jo dr JO ^ 



2n 



u' + 



n*AA 



2tt/uji 27r , 

= — u. 

C^i 



(84) 



UJi 

Here m* denotes the constant value (|76|) of dt/dr for the circular reference 
orbit. At this time the values of the angular co-ordinate of the perturbed and 
circular reference orbits coincide. The angular direction of the periastron can 
therefore be calculated by calculating the angular co-ordinate of the circular 
orbit at time T. The angular shift of the periastron per period is then 

Sip = ^{i = T) - ^(0) - 27r, (85) 

where we have used overlines to denote the co-ordinates of the circular orbit. 
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Now as u"^ is conserved, we have 



^(^) - ^(0) = 4§ = ^- (86) 



R^ w u 
Substitution of this resuh into eq. (PH) , using (RO) , gives 



This is in complete agreement with eg . (|6lD . 

5 Conclusion 

The method of investigation of particle motions close to the exact solutions 
already known, can be used in other exact solutions of Einstein-Maxwell 
system. We believe that it is worthwhile to pursue the efforts in this direction, 
in particular because this method introduces the small parameter (the norm 
of the deviation vector) different from the small parameters used by many 
authors in search for approximate solutions of the theory, such as e.g. the 
ratio M/ R or Q/ R. It is also well-suited for computer calculations. 
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A Appendix: Connections and curvatures for 
Reissner-Nordstr0ni geometry 

In this appendix we collect the expressions for the components of the con- 
nections and Riemann curvature used in the main body of the paper. 
a. Connections. From the line-element (^) one derives the following expres- 
sions for the connection co-efficients: 

,^ Mr-Q' 

'^ '' r (r2 - 2Mr + Q^) ' 



r^/ = ^ (Mr - Q2) (r^ - 2Mr + Q^ 



r,; = sin^ e r,/ = - ^ (r^ - 2Mr + Q^) , (88) 

r ^ = r "p = - 

re '- np ^1 

^e/ = —R^ r ^ = - sin ^ cos Q. 

^ sm d ^^ 

b. Curvature components. The corresponding curvature two-form compo- 
nents R^u = I R^x^ydx'^ A dx^ are: 

R^^ = — {2Mr - 3g2) dt A t/r, 

i?ie = - ^ (Mr - Q^) (r^ - 2Mr + Q^) dt A rf^, (89) 

Rt^ = -^ (Mr - Q^) (r^ - 2Mr + Q^) sin^ 6 dt A d<^, 
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and 

Mr - Q"^ 



«^.-;.^|^-^''*-''- '^°^ 



Rq^ = - i^2Mr - Q'j sin^ 9d9 A d^. 
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